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ABSTRACT 

Observations of spiral galaxies show a strong linear correlation between the 
ratio of molecular to atomic hydrogen surface density -Rmoi and midplane pres- 
sure. To explain this, we simulate three-dimensional, magnetized turbulence, 
including simplified treatments of non-equilibrium chemistry and the propaga- 
tion of dissociating radiation, to follow the formation of H2 from cold atomic gas. 
The formation time scale for H2 is sufficiently long that equilibrium is not reached 
within the 20-30 Myr lifetimes of molecular clouds. The equilibrium balance be- 
tween radiative dissociation and H2 formation on dust grains fails to predict the 
time-dependent molecular fractions we find. A simple, time-dependent model 
of H2 formation can reproduce the gross behavior, although turbulent density 
perturbations increase molecular fractions by a factor of few above it. In con- 
tradiction to equilibrium models, radiative dissociation of molecules plays little 
role in our model for diffuse radiation fields with strengths less than ten times 
that of the solar neighborhood, because of the effective self-shielding of H2. The 
observed correlation of -Rmoi with pressure corresponds to a correlation with local 
gas density if the effective temperature in the cold neutral medium of galactic 
disks is roughly constant. We indeed find such a correlation of -Rmoi with density. 
If we examine the value of -Rmoi in our local models after a free-fall time at their 
average density, as expected for models of molecular cloud formation by large- 
scale gravitational instability, our models reproduce the observed correlation over 
more than an order of magnitude range in density. 

Subject headings: astrochemistry — molecular processes — ISM: molecules — 
ISM: clouds 
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Introduction 



Stars are universally observed to form in molecular clouds. Indeed, recent observations 
have demonstrated that the surface density of star formation Sstjr co r relates linearly with 



the surface density of molecular hydrogen Sh ., (IRownd fc Young Il999l : IWong fc Blita 12002 



Gao fc Solomonll2004l : iBigiel et al.ll2008l . l201l[ ) . This has led to the suspic i on that the forma- 
tion of stars depends on the formation of molec ular hydrogen (jSchayell2004j : iKrumholz fc McKee 
20051 : lElmegreerul2007l : lKrumholz et al.ll2009bl ) . However, the opposite view has also been ar- 



guec 



1999 



2007 



that molecular clouds and then stars form from converging flows (IBallesteros-Paredes et al. 



Kovama fc Inutsu ka 2000';^ 



Hartmann et al. [20011 



Hennebelle et al.lboOSt iHeitsch k HartmamJIgOOsi 



drive r i by large-scale gravitational instability (iRafikov 



2002 



Kravtsov 



2003 



Shetty &: Ostrikei 



Dalcanton et al. 



2008 



2004; Mac Low fc Klessen 



Ostriker et al 



Vazguez-Semadeni et al. 



Inoue fc Inutsuka 



2001 



Kim fc Ostriker 



2004 



2006 



2009 ) , primarily 



200 



Li et al 



Ballesteros-Paredes 



Elnaegreen 



2005 



2006 



2OIOI ). In this picture, the formation of molecules is 



a consequence, not a cause, of the conditions required to form stars. 



Wong k Bhtzl (l2002f l and iBhtz k Rosolowskvl f l2004j . 120061 ) added to the puzzle with 
the observation that in the inner parts of spiral galaxies the ratio of molecular to atomic 
hydrogen surface density 

-Rmol = Shz/ShI OC Pm-, (1) 

where Pm is the midplane pressure and Shi is the surface density of H I. In detail, they 
derived an empirical pressure dependence for -Rmoi that can be fit se veral different ways. 
Takin g a common fit to all their individual data points yields the fit (iBlitz k Rosolowsky 
2006h 

Pm I k 



0.94±0.02 



[(4.5 ± 0.14) X lO^J ' 

where Pm/k has units of K cm~^. This result was extended by iLeroy et al.l (120081 ) who 
studied a broad sample of spiral and dwarf galaxies from THINGS, and gave the results in 
their Table 6. They found results quite consistent with equation ([2]). 



Both iBlitz k RosolowskyI ( l2006l ) and iKrumholz et al.l ( l2009al ) offered explanations of 
this observed correlation that rely on the equilibrium balan ce between radiative dissoc iation 
and H^ molecule formation on the surfaces of dust grains. iBlitz k Rosolowskvl (120061) built 



on the original suggestion by lElmegreenI ( 1l989l ) and lElme green k Parravand (119941 ) that the 
midplane pressure in disk galaxies de termines molecular cloud properties. They began with 
the calculation by lElme green! (jl993[ ) that the equilibrium fraction of gas in the molecular 
phase 

/mol = SH2/Stot (3) 
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depends on both the interstellar pressure and the radiation field j as 

/mol j 



(4) 



This equation was derived by assuming two populations of spherical clouds, one with con- 
stant density, representing diffuse clouds, and one w ith dens ity profiles, represent- 



i ng hydrostatic, self-gravitating c l ouds. Observations jHeilesI l200l[) and turbu lent models 
ke Avillez fc BreitschweTdtl boosl : Ijoung fc Mac Lowl booei Iciover et all boioh reveal that 
the distribution of densities among the diffuse clouds is continuous rather than the delta 
functions predicted by thermal phase transitions, demonstrating some of the limitations of 
this approach. 



For low values of -Rmoi ~ 0.5, the approximation Rmoi 



/mol can be made. iBlitz fc Rosolowsky 



( 120061 ) then invoked the observed correlation Ssfr oc , and assumed that the star forma- 
tion rate determines the local radiation field directly (j oc Ssfr) to derive that 



Pmo\ Pf 



1.2 



(5) 



However, the discussion in iBlitz &: RosolowskyI ( 120061 ) and iLeroy et al.l ( 120081 ) neglects that 
this relationship does not hold for the high values of -Rmoi S> 0.5 observed in some galaxies, 
as Equation (jl]) does not hold, nor does the approximation Rmoi — /moi- 



Krumholz et al.l (j2009al Jbl). and iMcKee fc Krumhola (120101 ) calculate the equilibrium 



value of /mol in a uniform- density, spherical gas cloud exposed to far ultraviolet (FUV) 
radiation. An es sentially identical result was already derived for a slab model of a cloud by 
Sternberg] ( 1l988l ). They all demonstrate that the fraction depends on only two parameters. 
The first parameter is the characteristic dust optical depth th = nncrdR of a cloud of radius 
R, where uh is the number density of hydrogen nuclei in the atomic gas, and ad is the 
cross section for dust absorption per hydrogen atom. The second parameter is the effective 
intensity of the ionizing radiation 



X = fdOdcEl/{nHTl), 



(6) 



(denoted aG by ISternberglll988[ ) . where fd is the fraction of absorbed photons in the Lyman- 
Werner bands that result in dissociation, cEq is the ambient flux of radiation in the Lyman- 
Werner bands, and TZ is the rate coefficient for the formation of H2 on grains. 

The effective intensity x can be written in terms of the value of Gq, the ratio of Eq 
to the typica l dissociating radiation fleld in the Milky way of Gq = 7.5 x 10~^ photons 
( lDraindll978l ). as 

X = 7lK_2l/7^_l6.5)(G';/r^^,), (7) 
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where ad-21 = o-rf/lO cm and 7^-i8. 5 = 7^/10 



-16.5 (,j^3 g i_ 



Both scaled constants are 



of ord er unity in the solar neighborhood. iKrumholz et al.l (l2009bl ) and iMcKee &: Krumholz 
(I2OIOI ) go on to express x terms of the metallicity and the physical constants only, by 
assuming that the atomic gas is in two-phase thermal equilibrium. However, this second 
step is not required. Not taking it allows application of their chemical model to a cloud of 
arbitrary density and ambient radiation field. 

These two papers then extend their approximation to compute /moi in a spherical cloud 
composed of both atomic and molecular gas, by making the assumption that the atomic gas 
is in pressure equilibrium with the colder, denser molecular gas. This allows them to derive 
a fit to /moi for the spherical gas cloud as a function solely of ^tot, the dust cross section 
ad, and the mean mass per hy drogen nucleus fin- Using the improved approximation given 
bv lMcKee fc Krumholzl J2010h . 



mol 



1.1 5 



Aj l + s/4' 



for s < 2, where, 

s = ln(l + 0.6x + 0.01x^)/(0.6r,), (9) 

and the optical depth through the atomic-molecular cloud Tc = 0.75'EtotO'd/ fJ'H- For s > 
2, fmni = is a bette r approximation. Essentially this same relationship was used by 



Krumholz et al.l ( l2009bl ) to derive the star formation rate in galaxies, taking into account 
the observed correlation between molecular hydrogen and star formation surface density. 

But why should H2 exert such a strong influence on star formation? Although H2 is a 
coolant, it is effective only down to temperatures of 200 K, whereas atomic fine structure 
emission (primarily from C"*") can cool the gas down to temperatures of around 60 K in 
regions of low dust extinction, or as low as 15 K in more highly shielded regions w here 
the photoelectric heating rate is small (iGlover fc Mac Lowll2007bt iGlover fc Clarkll201l[ ). To 
cool further, to the 10 K typical of most prestellar cores, generally requires CO, which indeed 
forms efficiently only in H2-rich gas. However, the difference between 20 K and 10 K can 
only slightly affect th e star formation rate, as demonstrated by the recent simulations of 
Glover fc Clm^ toilh . 



On the other hand, gr avitational instability produces dense gas that quickly forms H2 
( IGlover &: Mac Lowl l2007bl ). Thus, H2, and other molecules that form with it, such as 
CO, may pr imarily just ac t to tr ace dense gas that is already gravitationally unstable and 



collapsing. lOstriker et al.l ( 120101 ) present an analytic model for how the combination of 



dynamical and thermal equilibrium may lead to this situation, also yielding a prediction of 
the correlation between the molecular surface density and the midplane pressure based on the 
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linear relation derived between radiation field strength and pressure. However, this model 
relies on the empirical correlation Ssfr oc without offering a theoretical explanation for 
it. 

In this paper we examine the evolution of the molecular hydrogen to atomic hydrogen 
ratio -Rmoi as a function of time and examine whether that evolution can explain its observed 
correlation with pressure, and thus its relationship to the star formation rate. To do this, 
we have incorporated a simplified, non-equilibrium, chemical network and radiative transfer 
model into a three-dimensional simulation of supersonic (and super- Alfvenic) , magnetohy- 
drodynamical turbulence. This all ows us to follow the formation history of H2, st arting from 
turbulent, magnetized atomic gas (iGlover fc Mac Lowll2007al : iGlover et al.ll2010| ). 



Our models follow the interplay between the cold neutral medium (CNM) and molecular 
gas at different average densities. They have velocity dispersions of 5 km s~^, with resulting 
densi ties dispersed arou nd the mean value by more than four orders of magnitude (see Fig. 
7 of iGlover et al.l 120101 ). Although we do not explicitly model gravitational collapse, we do 
examine regions with varying average density, providing a first approximation to the sequence 
of quasi-static states that a large-scale, gravitationally contracting region will pass through. 
Only a small fraction of the molecular gas will ev er undergo full-scale gravitational collapse 
(iGoldreich &: Kwanlll974l : iKrumholz fc Tanll2007l ). with the rest remaining at pressures close 
to the surrounding atomic gas. Because we do include the full suite of heating and cooling 
processes relevant for both cold atomic and molecular gas, our models offer insight into the 
behavior of the cold gas even on large scales. 

At any particular time after molecules begin forming, we find that the molecular fraction 
-Rmoi varies primarily as a function of local gas volume density, almost independently of 
radiation field strength, except in diffuse regions with mean extinctions of Ay ^ 0.3. This is 
because the strong self-shielding of H2 means the rate-limiting step in denser regions is the 
slow formation of molecules on dust grains. Significant m ass fractions of the flow in these 
models become fully molecular in only a few million years (IGlover Sz Mac Lowll2007bl ). but 
molecule formation continues over periods exceeding 20 Myr, comparable to cloud lifetimes. 
Thus, equilibrium values of the molecular fraction cannot be relied on. This is contrary to 
the behavior of CO, whose fractional abundance is far more sensitive to the dust column 
density (and thus radiation field strength) than to th e volume density, because of its much 
weaker ability to self-shield (IGlover fc Mac Lowll201ll ). 



If molecular clouds rarely or never reach equilibrium, then what determines the observed 
molecular fraction? The hypothesis that global star formation is controlled by gravitational 
instability of the gaseous disk has been consid ered at least since the calc u lation of the grav- 
itational instability criterion in such disks by IGoldreich fc Lynden- Belli ( 119651 ). Numerical 
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simulations by iKravtsovl (|2003[ ): iLi et al.l (|2006[ ). and lTasker fc BryanI (120061 ) sup ported this 
idea, demonstra ting that it naturally could explain observati ons of both global ( iKennicutt 
19981 ) and local ( iMartin fc Kennicuttll200ll : iBigiel et al.l 120081 ) correlations between gas sur- 
face density and star formatio n rate. The clunapy nature of high-redshift galaxies c an also be 
explained by this hypothesis (IBournaud et al.ll2007l : iBournaud fc ElmegreenI 120091 ). Molecu- 
lar cloud formation as a by produ ct of gravitationa l instability has been argued to be a natural 



consequence of this hypothesis (lElmegreen 



2002 



Ballesteros-Paredes et al.ll2007l : iTasker fc Tan 
the free-fall time 



20091 



Mac Low fc KlessenI 12004 iLi et al.l 12005 
}). This suggests that the key timescale is 



-1/2 



(10) 



at the local average density p =< Hunu >■ We show that our non-equilibrium models 
indeed predict the observed molecular fractions to occur after a free-fall time at the average 
density in each model. 

In section |2] we describe our numerical method and the models that we ran. We then 
offer a simple analytical model of molecular hydrogen formation in section |3] and discuss how 
it succeeds and fails in comparison to the full numerical results. In section H] we compare 
our results to equilibrium models and the observations, while in section |S] we offer caveats 
and consider the implications of our results. 



Numerical Method 



To study the formation of molecular hydrogen from atomic hydrogen, we use mesoscale 
simulations of turbulent, magnetized, atomic gas incorporating a non-equilibrium chemical 
network, including an approximation to the background dissociating radiation field. In order 
to resolve the turbulent flow sufficiently well to follow the formation of H2, we neglect the 
large scale structure of the interstellar medium, and instead simulate periodic boxes with side 
length L = 5-20 pc having defined average initial number density no (and other parameters 
such as metallicity Z, turbulent velocity dispersion, and magnetic fiux). 



We use a version of ZEUS-MP (IHayes et al.ll2006l ) modified to include a subcycled chem- 
ical model and the ass ociated radia t ive an d chemical heating and cooling. The method was 
described in detail in I Glover et al.l ( I2OIOI ) but for the models presented here we explore a 
wider range of initial conditions and run several of the simulations for considerably longer. 
Radiative transfer is included usi ng a six-ray approximation , with column densities mea- 
sured from the edges of the cube (iGlover fc Mac Lowl l2007al ) . Although this is not strictly 
self-consistent, it does capture the basic dynam ics of radiative dissociation well (also see 
Glover et allboiol and blover &: Mac Lowlboilh . 
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All of the mo dels presented here i nclude driven turbulence, rather than the decaying 
turbulence used in iGlover fc Mac Lowl (l2007bl ). with an rms velocity Vrms = 5 km s~^, and 
an initial vertical magnetic field of Bq = 5.85 fiG. In Tabled] we give the initial density uq 
and metallicity Z, the box size L, the strength of the radiation field G in units of Gq, and 
the number of zones along one edge of the cube nx. 



Molecule Formation 



Our models show that the molecular to atomic ratio -Rmoi grows almost as a power 
law in time over many dynamical times (Figured]). Although our models will eventually 
reach saturation, the time taken to do so is generally longer than the plausible lifetimes 
for molecular clouds of 20-30 Myr (iFukui fc Kawamural l2010l ). We only find evidence for 
reaching an equilibrium value of -Rmoi within the first 25 Myr of evolution of our models for 
the model with the least active molecule formation (the low density and metallicity model 
n30-Z01 discussed in 5 0). 



3.1. Formation Time 

We can derive a simple analytic model for the behavior seen in our simulations by 
considering the finite formation times cale of molecu l ar hy drogen in a uniform medium. A 



rather similar approach is followed by iGnedin et al.l (120091 ) . The fractional local density of 
molecular hydrogen 

= 2niijn, (11) 

where n = nn + 2nH2 ? grows as 

Xn, = C{t)nun{T), (12) 

where the formation rate TZ oc T^/^, and we have included a time dependent clumping factor 
C to account for the density perturbations. We can eliminate nn from these two equations 
to find 

Xu, = 2Cnn,TZ{l/Xu, - 1). (13) 

Integrating this equation with a change of variable and the assumption of a uniform value 
for the clumping factor C, we find 



= 1 - exp{-t/TF), 



(14) 




t (Myr) 

Fig. 1. — Time history of the value of the column density ratio -Rmoi integrated through 
turbulent, magnetized boxes, starting with purely atomic gas, for models listed in Table [1] 
Models shown have average number density 30 (purple, lowest), 100 (blue), 300 (green), and 
1000 (brown, highest) cm~^. Canonical models (solid) are compared with GO models lacking 
radiation (dashed), and -L5 models with smaller boxes (dotted). The thin, solid, black curves 
show the analytic model given by Equation ( IT^ for time constants given by Equation (fT5l) 
for the densities of each model, while the thin, dashed, black curves show the analytic model 
with a clumping factor of two. The small box models diverge from the unclumped predictions 
earlier, as expected. No single clumping factor exactly reproduces all results. 
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where = 1/ [CTZn) is the time constant for H2 formation. iHoUenbach fc McKed (119791 ) 
aheady showed that for typical molecular cloud temperatures and without clumping (C = 1), 



1 Gyr 



n/{l cm 



(15) 



The assumption of no clumping should be good for times much less than the dynamical time 
tiiyn = L/vrms, whcu the average density rather than the peak clump density determines the 
evolution. In our canonical model, tdyn = 3.9 Myr, while our small box models with suffix 
L5 have tdyn = 0.98 Myr. 

We can compare this simple model with our numerical results and the observations 
under the assumption that the local density ratio X^^ can be used to derive the column 
density ratio i?moi = -'^H2/(1 ~ -^H2)- Figure [T] shows that, as expected, this analytic model 
works well at early times, up to roughly 0.4 Myr for the canonical model, and 0.1 Myr for 
our small box models. In both cases, this represents a clumping time 



6tr ~ O.lt 



dyn • 



(16) 



Thereafter, Rmoi increases above the predi ction of the simplest mode l, because of the en- 
hanced clumping caused by the turbulence (iGlover &: Mac Lowll2007bl ). 



3.2. Clumping 



The assumption of an equilibrium fraction of molecular hydrogen fails, not just because 
of the long formation timescale of molecular hydrogen but also because molecular hydro- 
gen does not form uniformly. Supersonic turbulence in the cold, neutral, atomic medium 
produces strong, transient, density perturbations. Peak densities exceed the average value 
by the square of the local thermal Mach number, which for conditions in the cold neutral 
medium can be a factor of 6n/n > 25. Since the formation timescale (Eq. [T^ is inversely 
proportional to density, formation proceeds far more quickly in the density peaks. Because 
of the strong self-shielding of H2, it does not immediately photodissociate when it subse- 
quently advects into lower density region s, leading to non-equilibrium fractions in lower 
density regions ( iGlover fc Mac Lowll2007bl ). 



No constant value of the clumping factor fits the late time results, however. For example, 
the value of C = 2 used in Figure [T] fits the n = 30 cm~^ results well at late time, but 
overpredicts the canonical n = 1000 cm~^ result substantially, though it does match the 
small box model well at late times. The time and densi t y depe ndent variation of the clumping 
factor has begun to be studied by iMilosavljevic et al.l (l201l[ ). 
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4. Comparisons 
4.1. Equilibrium Models 

We compare a selection of our tim e-dependent results fo r turbu lent gas to the equilibrium 



values derived for spherical clouds by iMcKee &: Krumholzl (120101 ) . In Figure [2] we compare 



time histories of a representative set of our models to the approximate equilibrium values 
of i?moi for the conditions of each of these models given by Equation ([H]), taking Etot to be 
the column density through the cubical simulation domain, and -Rmoi = /moi/(l — /moi) The 
equihbrium approximation works particularly poorly at Rmoi ^ 1, predicting no molecules 
in many cases where we find measurable molecular fractions. At larger values of -Rmoi, our 
models do intersect the spherical, equilibrium values, but at widely varied times, and with 
no indication of asymptotic approach to those values. 

The equilibrium models predict strong variation depending on the ambient radiation 
field, while we see little to no effect (comparing the GO and GIO models to the canonical 
models, we see that they are practically identical in the n300 case). This is because the 
evolution is controlled by H2 formation rather than photodissociation equilibrium. This 
occurs because the supersonic turbulent flows have strong density enhancements, leading 
both to efficient self-shielding in the regions where most H2 is forming and also fast local 
molecule growth in those enhancements. 

We do note that the equilibrium models in the intermediate regime give the qualitatively 
correct result of H2 formation at moderate densities. Although this is neither quantitatively 
corrrect nor extendible to other regimes, it nevertheless means that observational compar- 
isons that rely only on this qualitative result can be successfully made to the equilibrium 
models. 

The value of -Rmoi in our models depends on the metallicity, through its regulation of 
the dust density, and thus the H2 formation rate. However, only our very lowest density and 
metallicity model n30-Z01 approaches an equilibrium value during the 20 Myr that we ran our 
models. This suggests that real molecular clouds rarely find themselves in equilibrium with 



the d iffuse radiation field because their lifetimes are not long enough (IFukui fc Kawamura 



20101). Instead, they are probably dissociated and dispersed by intense, local, UV radiation 



from newly formed stars. 
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Fig. 2. — Time history of -Rmoi for models with average number density 30 {black, grey), 
and 300 (green, blue) cm~^. Canonical models (solid) are compared with GO models lacking 
radiation (dashed), GIO models with ten times higher radiation (dash- dotted), L5 models with 
smaller boxes (dotted), and low metaUicity models with 0.3 (blue dot-dashed) and 0.1 (blue, 
grey dashed) the canonical metaUicity. The models with varying levels of radiation in the 
high density case almost entirely overlap each other. Thin ho rizontal lines and arrows sh ow 
the predictions from the equilibrium calculation described by iMcKee fc Krumholzl ( 120101 ) as 
given in Equation ([8]). The predictions for the radiation-free GO models are -Rmoi — ^ oo, 
while the high-radiation, low-density well as a number of other models, is predicted 

to have -Rmoi = 0, as indicated by the arrows at the top and bottom of the plot. The very 
weak dependence of our models on radiation strength directly contradicts the prediction of 
the equilibrium models. 
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4.2. Observations 



Averaged over kiloparsec scales, the average temperature, magnetic pressure, and tur- 
bulent velocity dispersion of the cold, neutral med .ium vary little, with T cnm ~ 60-70K 
20031 ) ■ field strength B ~ 5/iG 



( iWolfire et al 



1995 



Hei 



higher as seen in observations of ex ternal galaxies (iBeckl 120051) , and turbulent Mach number 
Mj ~ 5 (IHeiles fc Trolandl l2003l ). lOstriker et al.l (120101 ) argue that this roughly constant 



es fc Trolandll2003[) or perhaps 20% 



effective dynamical temperature comes from self-regulation of star formation. The roughly 
constant dynamical temperature suggests that the apparent dependence of the fraction -Rmoi 
on pressure Pm given in Equation (HI) rn ay actually refiect a density dependence, as already 
suggested by iBlitz fc RosolowskyI ( 12004 ) . 



Figure [T] indeed shows a clear density dependence, but it also shows strong time depen- 
dence, raising the question of what molecular fraction should actually be compared with the 
observations. If the hypothesis holds that molecular clouds actually form during large-scale 
gravitational instability as discussed in § [H then the characteristic time in their lives at 
which clouds will be observed is a free-fall time (Eq. [TOj) after gravitational collapse begins 
in cooled atomic gas. Our models do not begin with clumped atomic gas as the actual ISM 
does. Instead, the turbulent fiows initialized on the grid clump the gas over a finite time 
6tc (Eq. [T^ . In order to account for this effect, we have taken the zero point in time for 
this comparison to be Stc rather than the initial time of our models. We emphasize that the 
observations show an order of magnitude scatter around the average RmoX^Pm correlation, 
so that using a characteristic time to describe a time-dependent process appears likely to 
be consistent with the observations if they catch clouds at different points in their evolution 
centered around a free-fall time. 

Each of our models has a fixed average density uh given in Table [H We use that value 
to compute the average pressure in our models 



P = nnkTcNM + (-B^)/87r + O.bfinnHvl 



(17) 



24 



where we assume that Tcnm = 60 K, and the mean mass per particle fin = 2.11 x 10 
appropriate for atomic gas with 10% helium by number. Because our models are subject to 
dynamo generation of magnetic field from the initial value, we use the actual values of the 
average magnetic pressure at each time rather than the initial value. 

We then use Equat i on (El) to derive the value and standard deviation of -Rmoi predicted 



by iBlitz fc RosolowskyI ( l2006l ) for the pressure for that mo del by setting Pm. = P- We also 



derive values of -Rmoi from the spiral and dwarf samples of iLeroy et al.l ( 120081 ). The spiral 
sample has Rmoi = (IQ-^-^^ Pm/kf-^°, while the dwarf sample has Rmoi = {lO'^-^^ Pm/kY'^^. 
All of these values are plotted in Figure |3] at times tfi{nH) after the beginning of evolution 
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of our models. 

We find that our models predict the observed -Rmoi across more than an order of mag- 
nitude in pressure, with an accuracy better than the scatter in the observationally derived 
relations themselves. The ability to predict not just a single value, but a whole family of 
values accurately suggests that we have captured the important physics of the variation of 
the molecular fraction _Rjnoi with the estimated midplane pressure. 



5. Discussion 

5.1. Influence of numerical resolution 

Figure [3] shows that the values of Rmoi we find in our L = 5 pc simulations are sys- 
tematically higher at late times than those found in our L = 20 pc simulations. (This is in 
addition to the early time differences discussed above that occur because of the difference 
in crossing time tdyn between the two box sizes.) Since both sets of runs were performed 
with the same numerical resolution of nx = 128 zones per side, the spatial resolution of our 
L5 runs is four times larger than that in our other runs. This prompts the concern that 
our simulations may be under-resolved, and that our reported results for -Rmoi niight have a 
significant dependence on the spatial resolution of the simulations. 

To investigate this, we have performed several additional simulations with a higher nu- 
merical resolution. In Figure Hl^a), we compare the evolution of -Rmoi in two canonical runs 
and the corresponding L5 runs, with that in two runs performed using the same sets of 
simulation parameters, but with twice the numerical resolution. Increasing the resolution to 
nx = 256 zones per side leads to slightly faster production of H2 at late times, particularly in 
the lower density runs, since we are better able to resolve the highest density structures pro- 
duced by the turbulence. However, the effect is small, and does not explain the significantly 
larger difference that we see between the results of the 5 pc (L5) and 20 pc runs. 

Ideally, we would perform a similar comparison using results from runs with L = 20 pc 
and a numerical resolution of 512^ zones, which would have the same spatial resolution as our 
L = 5 pc runs. Unfortunate l y, the high computational cost of performing a 512^ simulation 



using the full iGlover et al.l (120101 ) cooling and chemistry model has so far rendered this 
impractical. 

However, it has been possible to perform simu lations with this resolut ion by adopting 



a simpler treatment of the CO chemistry based on iNelson fc Langerl (119991 ). By using this 



simpler treatment, we accept a slightly higher degree of error in our predicted CO, C and 




t (Myr) 

Fig. 3. — Time history of the value of the column density ratio -Rmoi for the same models 
presented in Figure [T], shown on semi- log ax es for clarity. Superp o sed, w e show the empirical 
fit for the canonical model of each density of iBlitz fc Rosolowskyi (120061), includii ig the errors 



fx), a s well as the spiral (diamond) and dwarf (triangle) galaxy samples of iLeroy et al. 
( 120081 ). For all points we take a CNM temperature of 60 K, and the average density, time- 



dependent magnetic pressure, and turbulent pressure of each model. They are placed on the 
figure at a free-fall time for the average density of that model, as described in the text. To 
account for clumping of the atomic ISM, we use a zero time of 6tc in this plot, so that the 
initial time of each model is negative. 
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abundances, and lose our ability to follow the evolution of chemical species such as OH 
or water, but still retain the ability to follow the evolution of the H2 abundance accurately, 
since our treatment of the hydrogen chemistry remains the same. At these densities, the 
change in cooling from the modified chemistry will produce negligible changes in the chemical 
evolution. We show the results of this full resolution study in Figure Hl^b) . Once again, we 
see that although our results remain somewhat sensitive to the numerical resolution of the 
simulation, the effect is relatively small and cannot be responsible for the difference between 
the L = 5 pc and L = 20 pc runs. 

If spatial resolution is not the answer, then what is responsible for the difference between 
the runs with different box sizes at late times? The lower visual extinction of the gas in the 
smaller box allows stronger photoelectric heating, which leads to the gas in the L5 runs 
having a systematically higher temperature than the gas in the larger volume simulations. 
For example, the volume weighted mean temperature of the gas at the end of run n300 is 
approximately 40 K, while in run n300-L5 it is 66 K. This difference in temperatures leads 
to a systematic difference in the H2 formation rate coefficient, owing to the T^/^ temperature 
dependence of the rate coefficient, and hence leads to a faster production of H2 in the smaller 
volume simulations. 



5.2. Caveats 

In this section, we discuss the limitations of this work. The biggest problem is that we 
are substituting a study of local isotropic turbulent regions without explicit self-gravity for 
one of the dynamical interstellar medium. However, to reach the resolutions of less than 0.2 
pc used here in any large-scale model of the interstellar medium remains computationally 
challenging. These resolutions are required in order to resolve the turbulent density peaks 
that determine the H2 formation rate. 

We approximate the external field of FUV dissociating radiation by injecting radiation 
at the boundaries of our periodic box. Although this is physically inconsistent, we find 
that in this study it makes virtually no difference, as models with and without FUV behave 
almost identically. This is because the densities that we study are sufficiently high that 
H2 self-shielding efficiently blocks penetration of UV into the dense regions where H2 forms 
and remains. This also means that the opening of voids in the density distribution by the 
turbulence has much less effect on the H2 fraction than the creation of density peaks. 

In our models, we do not include local sources of dissociating or ionizing radiation. Star 
formation of course produces such local sources, which can be effective at dispersing molec- 
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ular clouds on short time scales. Turbulent models of such d ispersal have n ot ye t been per- 



forme d, but analytic models assuming spherical symmetry by iMatznerl (120021 ) and lElmegreen 



( 120071 ) show that large clouds can be dispersed in a dynamical time. Similarly, analytic mod- 
els including accretion sho w lifetimes comparable to the observed values of around 20 Myr 
(iFukui fc Kawamurall2010l ) in e ither a slab converging flow (jZamora-Aviles fc Vazquez-Semadeni 
201ll ) or a spherical geometry (iGoldbaum et al.ll201ll ). These relatively short cloud destruc- 
tion time scales support our assumption that most clouds are observed at an age of roughly 
a free-fall time. 

Our chemical model has been simplified to allow computation on a three -dimensional 
grid. However, in the range of temperatures and densities considered here, iGlover et al. 



(I2OIOI ) have shown that we reproduce the results of the UMIST model to within a few tens 
of percent. Moreover, we do not expect any of the chemical processes omitted from the 
model to have a significant impact on the fractional abundance of molecular hydrogen, and 
hence our values for Rmoi should be quite accurate. 

Finally, we assume in this study that the dominant contributions to the observed column 
densities of atomic and molecular hydrogen come from cold, dense clouds, rather than from 
the warm, diffuse medium surrounding these clouds. This is good to better than a factor 
of two in the Milky Way where cold gas rep resents rough ly 6Mq pc~^ while warm atomic 
and ionized gas is a bit less than 5Mq pc~^ (jFerrierell200ll ). and is probably better in more 
massive disks, but will break down in outer disks or other regions where the WNM component 
dominates the mass budget. Since in practice Rmoi is only measured in regions where CO 
emission can be detected, our approximation should be reasonable in the regime where 
observations are available. 



5.3. Comparison to Other Tests 



Krumholz &: Gnediru (1201 if ) tested the equilibrium models of iKrumholz et al.l (l2009b( ) 
(see eq. IH} against relatively low-resolution (smallest zone size Ax = 65 p c) global numerical 



simula tions of H2 formation in galaxies. T he simulations ^ yere d e scribed by iGnedin fc Kravtsov 
( I2OIOI ). and used the H2 formation law of iGnedin et al.l ( l2009l ). iKrumholz fc Gnediru (j201l[ ) 
found good agreement in cases where we find that our models reach high molecular fractions 
quickly (those with high density and low radiation intensity), but poor agreement in c ases 
where we find that time-dependence matters. The formation law of IGnedin et al.l ( 120091 ) ef- 
fectively reproduces the simple model that we have described in Equation (JH]). The inability 
of global models to resolve the density peaks is captured in the higher value of the clumping 
factor they require to reproduce observations. Because lower mean densities are deduced. 
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these models do still overestimate the sensitivity of molecule formation to radiation. 



Fumagalli et al.l (|2010[ ) compared the equilibrium model given by Equation to local 
(100 pc beam size) observations of the column density of H I and global (1 kpc beam 
size) observations of CO and H I. They demonstrate that the equilibrium model with freely 
adjusted clumping factor predicts high molecular fractions in regions of high local density, as 
we also predict without such an adjustable factor, and that there is a metallicity dependence, 
again as we also predict. 



5.4. Implications for the Star Formation Rate in Galaxies 

We have demonstrated that the molecular fraction -Rmoi remains strongly time dependent 
over periods of rather more than 20 Myr (Fig. [1]). Our high- resoluti on numerical results 



provide calibration for the time-dependent molecule formation model of iGnedin et al.l (120091 ) 
which assumes similar physics, but must rely on much larger clumping factors to reproduce 
sub-grid scale turbulence. Equilibrium values of -Rrnni are only reached at times well beyond 



the lifetimes of real molecular clouds of 20-30 Myr (iFukui fc Kawamurall2010l ). Therefore, 
explanations of the observed correlation between midplane pressure and -Rmoi that rely on 
equilibrium values seem unlikely to be correct. 

We find that the observed correlation can nevertheless be explained if most molecular 
gas occurs in regions that have been forming H2 for a free-fall time at their current number 
density (Fig. [3]). This is a natural prediction for the gravitational instability model for 
star formation in galaxies. Molecular cloud formation occurs as part of the same process of 
large scale collapse that eventually forms stars, so the observed correlation between molecular 
hydrogen and star formation rate surface densities comes about because they have a common 
cause, not because the first controls the second. 

Our result appears almost independent of the strength of the far ultraviolet dissociating 
radiation field during formation of the molecular clouds, as can be seen by the close agree- 
ment of our canonical models with otherwise identical models run without or with stronger 
dissociating radiation. In future work the upper limits of this behavior must be examined, 
since clearly a strong enough radiation field will eventually dissociate molecular gas, or at 
least ionize enough of it to dynamically disperse it, allowing dissociation to occur. In par- 
ticular, once massive star formation begins, the intense local radiation field probably does 
dissociate or disperse the remaining molecular hydrogen, terminating molecule formation 
relatively quickly and determining the cloud lifetime. 



The metallicity clearly influences the molecular hydrogen formation rate in our mod- 
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els, through its regulation of the dust density (Fig. [2]). The Hg formation rate depend s on 
the product nZ of the number density and the metallicity (iGlover fc Mac Lowl 1201 ll ) be- 
cause of the role of dust grains in the formation reaction, so this is expected. However, 
we find a weaker sensitivity than suggested by Equation derived from an equilibrium 
model. Detailed comparison of our results to galaxies with varying metallicity should prove 
informative. 
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Fig. 4. — Resolution studies of the time history of -Rmoi- (a) Comparison of runs nlOO (black, 
solid line) and n300 (blue, solid line) with runs nlOO-256 and n300-256 (dashed lines) having 
twice the numerical resolution, and with runs nl00-L5 and n300-L5 (dotted lines), having 
four times the effective resolution but in a four times smaller box. (b) Resolution studies 
of the nlOO (black) and n300 (blue) models reaching 512^ resolution by using the reduced 
chemical network of iNelson fc Langerl (119991 ). Runs ending in NL-128 (solid lines), NL-256 
(dashed lines), and NL-512 (dash-dotted lines) are directly compared to NL-L5 runs (dotted 
lines) that have physical resolution equivalent to the NL-512 runs. 
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Table 1. Input parameters used for each simulation. 



Run 


no 


z 


L 


G 


nx 




(cm j 


(^0) 


(pc) 


(<^o) 


(zones) 


n30 


30 


1 


20 


1 


128 


n30-G0 


30 


1 


20 





128 


n30-G10 


30 


1 


20 


10 


128 


n30-L5 


30 


1 


5 


1 


128 


n30-Z01 


30 


0.1 


20 


1 


128 


nlOO 


100 


1 


20 


1 


128 


nlOO-GO 


100 


1 


20 





128 


nl00-L5 


100 
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128 


nlOO-ZOl 


100 


0.1 


20 
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128 


nlOO-256 


100 


1 


20 
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256 


nlOO-NL-128 


100 
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20 


1 


128 


nlOO-NL-256 


100 
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20 
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256 


nlOO-NL-512 


100 
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20 


1 


512 
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1 


5 


1 


128 


n300 


300 
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20 
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128 


n300-G0 


300 


1 
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n300-G10 


300 
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20 


10 


128 


nQnn T K 
noUU-LiO 


oUU 


1 
i 





1 
i 


izo 


n300-Z01 


300 


0.1 


20 


1 


128 


n300-Z03 


300 
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20 


1 


128 


n300-256 


300 
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20 
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256 


n300-NL-128 


300 
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20 


1 


128 


n300-NL-256 


300 


1 


20 


1 


256 


n300-NL-512 


300 


1 


20 


1 


512 


n300-NL-L5 


300 


1 


5 


1 


128 


nlOOO 


1000 


1 


20 


1 


128 


nlOOO-GO 


1000 


1 


20 





128 


nl000-L5 


1000 


1 


5 


1 


128 
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Note. 



Zq is the solar metallicity, while Gq is the 



str ength of the l ocal interstellar radiation field as derived 
by iDraind ( 119781 ) . NL models were run without full carbon 
chemistry (see section [5?T|) . 



